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We utilize a fractional exclusion statistics of Haldane and Wu hypothesis to study the ther- 
modynamics of a unitary Fermi gas trapped in a harmonic oscillator potential at ultra-low finite 
temperature. The entropy per particle as a function of the energy per particle and energy per par- 
ticle versus rescaled temperature are numerically compared with the experimental data. The study 
shows that, except the chemical potential behavior, there exists a reasonable consistency between 
the experimental measurement and theoretical attempt for the entropy and energy per particle. In 
the fractional exclusion statistics formalism, the behavior of the isochore heat capacity for a trapped 
unitary Fermi gas is also analyzed. 
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I. INTRODUCTION 



The thermodynamics of a strongly interacting fermion 
system attracts much attention, either experimentally 
or theoretically, in the many-body community in recent 
years For the dilute contact interacting atomic sys- 
tem, the scattering length a describes the mutual interac- 
tion strength between two fermions. With the Feshbach 
resonance technology, it is possible to tune the s-wave 
scattering length a through changing the external mag- 
netic field. The scattering length a can be changed from 
a negative value to the other side of the resonance where 
a becomes positive. In the weak attraction Bardeen- 
Cooper-Schricffcr (BCS) regime with kp\a\ <§; l(with the 
Fermi wave vector fcp ), the mean field theory can give a 
reasonable description. With the increase of the interac- 
tion strength, the composite bosons can be formed in the 
fermion system, where the Bose-Einstein condensation 
(BEC) can occur [IS]. 

In the regime from negative scattering length to posi- 
tive value, the two-body interaction strength can be sin- 
gular with the existence of a zero-energy bound state. 
The limit kp \a\ — !• oo is called BCS-BEC crossover, which 
is also called the unitary regime The corresponding 
gas in this limit is called unitary Fermi gas |2|. At uni- 
tarity, the gas both has some properties of fermions and 
some of bosons, that is to say, the behavior of this gas 
is between that of fermions and bosons. The fermion 
system will show a universal thermodynamic behavior 
near the resonance point where the scattering length a 
diverges [H, 0]. 

Initially, physicists concentrated on the energy density 
per particle of the unitary Fermi gas at zero temperature 
theoretically. For a homogeneous gas, the ground-state 



energy per particle is given by E/N = liE F with the 
ideal Fermi energy Ep and the universal coefficient £ ■ 
It is E/N — \\f£,Ep for a trapped unitary gas, where 
Ep is the Fermi energy of the trapped noninteracting 
Fermi gas Determining the dimensionless constant £ 
has been an important physical subject in the past few 
years. 

Furthermore, exploring the thermodynamics of a uni- 
tary Fermi gas at finite temperature in theory to explain 
the experimental results is a more challenging task. A 
diagrammatic determinant Monte Carlo method for the 
negative-?/ Hubbard model was used to calculate the 
finite-temperature thermodynamics of a homogeneous 
unitary gas (4| . Through the quantum Monte Carlo sim- 
ulation, the finite-temperature thermodynamics of a ho- 
mogeneous unitary gas can be given by Refs.Jl, @. A 
self-consistent theory based on the combined Luttinger- 
Ward-De Dominicis-Martin variational formalism was 
also used to calculate it for a homogeneous unitary gas 

0- 

For the finite-temperature thermodynamics of a 
trapped unitary Fermi gas, the entropy and energy of 
the trapped unitary Fermi system have been experi- 
mentally measured in Refs.[§-fl0|. The quantum Monte 
Carlo simulation can also be used to calculate the finite- 
temperature thermodynamics of a trapped unitary gas 



11]. Various theoretical attempts have been established 
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to give the calculations for a trapped unitary Fermi gas. 
The pseudogap theory is one of them [T3, El ■ A mean- 
field method by solving the Bogoliubov-de Gennes equa- 
tions with an efficient and accurate method was used 
in Ref. 14]. A T-matrix calculation with a modified 
Nozieres and Schmitt-Rink approximation was adopted 
to explore it [lj|. The combined Luttinger-Ward-De 
Dominicis-Martin variational formalism was also proved 
to be a valuable approach for it [l6j]. We note that 
there are still differences among in these attempts on the 
finite-temperature thermodynamic properties of a uni- 
tary Fermi gas. 

Due to the scale invariance at unitarity, the thermo- 
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dynamic properties of the universal strongly interact- 
ing unitary fermions are related to those of the non- 
interacting ideal fermions. This means that the dynam- 
ical details will not appear in the final thermodynamic 
analytical expressions. From the microscopic point of 
view, the unitary fermions system is in between the 
fermionic and composite bosonic phases. To character- 
ize this crossover or intermediate unitary fermionic ther- 
modynamics, Haldane-Wu fractional exclusion statistics 
was used to discuss the finite-temperature unitary Fermi 
gas thermodynamics [13, [Hi- In physics, the thermody- 
namic behavior of the fractional exclusion statistics is be- 
tween those of Bose-Einstein and Fermi-Dirac statistics 
[T9I I20I ] , which is quite similar to the thermodynamics of 
real fermions at unitarity. As a hypothesis, the strongly 
interacting unitary gas is modeled by the non-interacting 
anyons. The priority is that the finite-temperature ther- 
modynamic properties can be investigated analytically 
plj . The quintessence hidden in this attempt is that the 
anyonic statistical parameter g characterizes the strongly 
interacting universal properties. 

Following the fractional exclusion statistics formalism, 
the aim of this work is to give a detailed discussion on the 
trapped gas thermodynamics. By generalizing the discus- 
sion for the homogeneous unitary system [22[ , we want to 
give the expressions and the corresponding numerical re- 
sults of entropy, physical chemical potential and isochore 
heat capacity of a trapped unitary system. The concrete 
comparisons with experimental measurements are made. 

The natural units fcs = h = 1 are used throughout the 
paper. 



II. FINITE-TEMPERATURE 
THERMODYNAMICS WITHIN FRACTIONAL 
EXCLUSION STATISTICS 

A. Framework of fractional exclusion statistics 

In this subsection, we will briefly review the formal- 
ism of the fractional exclusion statistics. By generalizing 
the simple formula with the fractional exclusion statis- 
tics, one can get the microscopic quantum states W of N 
identical particles occupying a group of G states [2(| HI] 



miGt-gNt-il-gW 



(1) 



where g is a statistical parameter, which denotes the 
number of states that one particle can occupy. For 
bosons, g — and g = 1 for fermions. 

If Gi and Ni are very large, one has the approximate 
expression of the logarithm of W 

In W ~ Y, K G i + (1 - 9)Ni) In (G. t + (1 - g)Ni) 

i 

- (Gi - gNi) In (d - gN) - N In N z ] , (2) 
through the Stirling formula In TV! = AT (In N — 1). 



The variational formulation of In W is 

d\nW =^[(l- ff )ln(G i + (l-< 7 )/V i ) 

i 

+ g]n(G i -gN i )-]iiN i \6N i . (3) 

However, SNi is not arbitrary, it must satisfy the con- 
ditions below: 



SN = J2 SNl = °> 

i 

5E = Y e i 6N i = °- 



(4) 



Setting the two Lagrange multipliers a = —ft/T and 
P = 1/T, one can have 

5lnW - aSN - PSE 
= ^[(l-^ln^ + (l-g)Ni) 

i 

+g In (Gi - gNi) - In N t + (ji - a)/T} SN t 
Gi 



E 



(1 - g) In ( 1 

Gi 



Ni ' 

(j, — e 
T 



SNi. 



(5) 



where [i is the chemical potential, T is the system tem- 
perature and ei is the single-particle energy for the state 
of species i. 

Defining the average occupation number Ni = N/Gi, 
through the Lagrange multiplier method 8 In W — aSN — 
/35E = 0, the most probable distribution of iV, can be 
derived from Eq.© 



Ni 



1 



(6) 



w< +g 

where u)i obeys the relation 

li-€i = -T [(1 - g) In (1 + wj) + gbxun] . (7) 

In Eq.([7]), we define uiq as the value of u)% at e% = 0. 
Consequently, one has 



fi = -T[(l -#)ln(l + uj ) + 3 In wo] 
At zero temperature, there is 

N = , for Ei > fx; 

N = - , for e, < n. 
9 



(8) 



(9) 



Furthermore, by inserting Eq.([B]) into Eq.([2]), one ob- 
tains the expression for entropy as 



S = lnW = Y J 



Gi 



w» + g 



[(u>j + 1) In (uji + 1) - uji hxuiip.0) 
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B. Thermodynamics of a trapped unitary Fermi 
gas 

The finite-temperature thermodynamics can be de- 
rived for the unitary Fermi gas trapped in a harmonic 
oscillator mm 2 r 2 /2, where the oscillator parameter w 
is defined as vj = (bJxLOyUiz) 1 / 3 and r denotes the po- 
sition of particles. The corresponding density of states is 
D(e) — e 2 jvo 3 . At zero temperature, the particle number 
and system energy are 



N = — 
9 



E = - 
9 Jo 



D(e)de 



eD{e)de 



El 
3m 3 '' 



4m 3 : 



(11) 



(12) 



where Ep obeys the relation Ep = g l / 3 Ep with the ideal 
Fermi energy Ep — (SN) 1 / 3 -^ in a harmonic oscillator. 

Comparing Eq. ([PTj) with Eq.(fT2"j) to eliminate tu, it is 
shown that 



E 



NEp 



-5 1/3 - 



(13) 



If the statistical parameter g in the fractional ex- 
clusion statistics is fixed through the zero-temperature 
ground-state energy, the finite-temperature thermody- 
namic quantities for a trapped unitary Fermi gas can 
be calculated. For the trapped system, the ground-state 
energy is related with the universal coefficient according 
to E/(\NE F ) = y/l d. From Eq.([13|), it is found that 
E/(^NEp) — g 1 / 3 . So the statistical parameter could 
be calculated as g — £ 3 / 2 = ^, with £ = | given by the 
developed quasi-linear approximation method |23T - [28j . 

For the finite-temperature trapped unitary Fermi sys- 
tem, the particle number and energy can be represented 
by turning the sum of quantum state into integral and 
changing the variable from de to dto 

n =5>#= r^=f^ 3 / l2 M, d4) 

; Jo u + g \wj 



E =y^GiNiei = 

Jo 



eD{e)de (T 



uj + g \vj 



where 



d/jj 



w(l + w) 



111 



1 + U) 

1 + UJQ 



1-9 



By replacing Eq. ([TT)l into Eq.([T4|) and Eq.(fl5|>. one gets 



3 (^) ft2 ^ = 1 ' 

E ( T \ 4 

_ = 3 (— J h 3 (uo), 



(16) 
(17) 



with the Fermi characteristic temperature Tp. 

One can turn the sum of Eq. pTj]) to an integral. With 
Eq. ([T6|) . one can get the explicit integral expression of 
entropy per particle 



S_ 

N 




(18) 



In order to derivate the expression of the isochore heat 
capacity, let us further discuss the particle number. From 
Eq. (TT4|) . the partial derivative of the particle number N 
to temperature T for fixed N and w is given by 



ON 



3T 2 

— 0-/12(^0) 



N,w 



dh 2 
dT 



dh 2 



dT 



N. 



From Eq. (fT9| . one can get 

d/j,\ fj, 3/i 2 (w ) 



dT 



N.zu 



T 2h 1 (uj )' 



(19) 



(20) 



Furthermore, with Eq. (|15j) and the thermodynamic re- 
lation between the isochore heat capacity Cy and internal 
energy E, one has 



dE 
df 



4| -) h 3 (cj )+[- 
vo \w 



T 



dh 3 
~df~ 



df 



(21) 



Substituting Ea.(|T4jl and Eq.® into Eq.([21]), one ob- 
tains the isochore heat capacity per particle 



Cy _ 4h 3 (oj Q ) 
N 



9h 2 {u] ) 



h 2 (u>o) 2/ii(w )' 



(22) 



III. DISCUSSIONS 

Based on the above analytical expressions, we will give 
the numerical results. 

The energy per particle versus the rescaled tempera- 
ture can be calculated from Eq. (fT6| and Eq. (fT7)) . As in- 
dicated by the Fig|TJ the energy for the trapped unitary 
Fermi gas with the statistical parameter <? = «f? versus 
the rescaled temperature is consistent with the experi- 
mental data [8J . 

Through Eq. fT?) and Eq. (fT5|) , the relation between the 
entropy per particle and the energy per particle is given 
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FIG. 1: The energy per particle versus the rescaled temper- 
ature. The solid curve denotes that for the trapped unitary 
Fermi gas, and the dashed one is that for the trapped ideal 
Fermi gas. The dots with error bars are the experimental data 
of Refs.fl. 
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FIG. 2: The entropy per particle plotted as a function of the 
rescaled energy per particle. The line styles are similar to 
FigH] The dots with error bars are the experimental data of 
Refs7[l,[icj. 
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FIG. 3: Similar to Fig[T] for the variation of the chemical 
potential as a function of the rescaled energy per particle. The 
dots with error bars are the experimental data of Ref . [Io| . 



imental data is not obvious. However, the chemical po- 
tential differs explicitly from the experimental result for 
E/(NEp) > 1.2. As the temperature increases, the 
experimental data also diverge from the curve of the 
trapped ideal gas. We can also plot the chemical po- 
tential versus the rescaled temperature in Figj4] 
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explicitly in Figj2j The entropy increases with increas- 
ing energy, and in the Boltzmann regime, the curve for 
a trapped unitary Fermi gas gets closer to and almost 
overlaps with that of the trapped ideal Fermi gas. It is 
found that the theoretical result is reasonably consistent 
with that of the experiment. 

The agreements in Fig[T] and Fig(5] between the theo- 
retical curves and the experimental data show that the 
statistical parameter g can be capable of describing the 
strong interaction of unitary Fermi gas at extremely low 
temperature. 

With Eqs.(0), JTBJ and ([17]), we have also shown the 
plot of the chemical potential varying with the rescaled 
energy per particle in Fig|3l It is a monotonically de- 
creasing function of the energy per particle. In the Boltz- 
mann regime, the curve for a trapped unitary Fermi gas 
gets closer to that of the trapped ideal one. In the ex- 
tremely low-temperature regime (E/(NEp) < 1.2), the 
departure between the theoretical result and the exper- 



FIG. 4: The relationship between chemical potential and the 
recaled temperature. The line styles are the same as in FigfT] 

From Eqs. ([TB)) and (|Tg)) , the entropy per particle as a 
function of the rescaled temperature is indicated in FigJS] 
In Fig|6l the isochore heat capacity per particle related 
with the rescaled temperature is presented. The numer- 
ical results are obtained by solving the coupled Eqs. p^)) 
and (|22l . 

Near the Boltzmann regime, the energy, chemical po- 
tential, entropy and isochore heat capacity of a trapped 
unitary Fermi gas are getting closer to and almost over- 
lap with those of the trapped ideal Fermi gas. In the 
low-temperature strong degenerate regime, the thermo- 
dynamic quantities of the unitary gas are lower than the 
ones of the ideal gas at the same temperature. The en- 
ergy, chemical potential and entropy of a unitary gas are 
always lower than the ideal ones at the same temperature. 
However, as shown in FigjBl the isochore heat capacity 
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FIG. 5: Same as in Fig|4]for entropy per particle plotted as 
a function of the recaled temperature. 
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FIG. 6: Plot of the isochore heat capacity per particle as a 
function of the recaled temperature. The line styles are the 
same as those shown in Fig[S] 



per particle given by the fractional exclusion statistics is 
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not a naive monotonously increasing function with the 
increase of the scaled temperature. This behavior is dif- 
ferent from that of the trapped ideal Fermi gas. 



IV. CONCLUSION 

The finite-temperature thermodynamic quantities of a 
trapped unitary Fermi gas have been discussed in terms 
of the approved fractional exclusion statistics framework. 
The study shows that the thermodynamic quantities of 
a trapped unitary Fermi gas will overlap with those of 
the trapped ideal Fermi gas in the Boltzmann regime. 
The thermodynamic quantities given by this formalism 
are lower than the trapped ideal ones at the same tem- 
perature except for the isochore heat capacity. 

The energy and entropy per particle manifest the con- 
sistency with the low-temperature experimental measure- 
ment. However, the detailed numerical study of the 
chemical potential shows that there is an explicit dif- 
ference between the result given in terms of fractional 
exclusion statistics and experimental data in the weakly 
degenerate regime. For the three-dimensional trapped 
unitary Fermi gas, the high-temperature weakly degen- 
erate behavior of chemical potential given by a simple 
fractional exclusion statistics hypothesis needs clarifying 
further. 
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